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Abstract. Inhomogeneous universe models have been proposed as an alternative 
explanation for the apparent acceleration of the cosmic expansion that does not require 
dark energy. In the simplest class of inhomogeneous models, we live within a large, 
spherically symmetric void. Several studies have shown that such a model can be 
made consistent with many observations, in particular the redshift-luminosity distance 
relation for type la supernovae, provided that the void is of Gpc size and that we live 
close to the center. Such a scenario challenges the Copernican principle that we do 
not occupy a special place in the universe. We use the first-year Sloan Digital Sky 
Survey-II supernova search data set as well as the Constitution supernova data set 
to put constraints on the observer position in void models, using the fact that off- 
center observers will observe an anisotropic universe. We first show that a spherically 
symmetric void can give good fits to the supernova data for an on-center observer, 
but that the two data sets prefer very different voids. We then continue to show that 
the observer can be displaced at least fifteen percent of the void scale radius from the 
center and still give an acceptable fit to the supernova data. When combined with 
the observed dipole anisotropy of the cosmic microwave background however, we find 
that the data compells the observer to be located within about one percent of the void 
scale radius. Based on these results, we conclude that considerable fine-tuning of our 
position within the void is needed to fit the supernova data, strongly disfavouring the 
model from a Copernican principle point of view. 
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1. Introduction 

The discovery of the dimming of distant type la supernovae (SNe la) [U |2] constituted 
the first impactful evidence that the expansion of the universe is in a phase of 
acceleration. This picture has been corroborated by several independent probes, 
including measurements of the cosmic microwave background (CMB) anisotropies [3] 
and baryon acoustic oscillations (BAO) |11[5]. Under the assumption that the universe 
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is homogeneous and isotropic, the apparent late-time acceleration is usually attributed 
to a mysterious energy component with negative pressure - dark energy - the nature of 
which is still unknown. 

In recent years, inhomogeneous universe models of varying degrees of complexity 
have been the focus of much attention as an alternative explanation of the apparent 
acceleration that doesn't require dark energy. The common base in these models is that 
the assumption of homogeneity and isotropy of the universe is an oversimplification and 
that it is the presence of inhomogeneities in the distribution of matter that gives rise to 
the apparent acceleration. Examples of inhomogeneous models include the Swiss-cheese 
models, which consider light propagation through a universe full of empty regions, and 
backreaction models, where cosmological perturbation terms are included in the overall 
dynamics of the universe (see [6] for a review of inhomogeneous models). 

In the simplest class of inhomogeneous models, we live within a large, spherically 
symmetric local void described by the Lemaitre-Tolman-Bondi (LTB) metric [Tf [8| [Qf [TO] . 
The LTB model is an exact solution of Einstein's equations containing dust only. In 
contrast to homogeneous universe models, the Hubble expansion rate and the matter 
density parameter depend on both time and the radial coordinate in LTB models. The 
apparent acceleration is caused by the radial dependence, such that our local underdense 
region has a larger Hubble parameter than the surrounding homogeneous universe. 

The LTB models violate the Copernican principle by placing the observer in a 
special place in the universe. Due to the observed near-perfect isotropy of the CMB, 
the observer is generally assumed to be located very close to the center of the void. 
Another challenge for the model is that the void must be of Gpc size in order to fit 
the SN la data, but such large voids are extremely improbable in standard models of 
structure formation [llj. The LTB models thus appear to be unrealistic, but based on 
observations they have not yet been ruled out. Several studies have shown that with 
an appropriately chosen void profile, the LTB models can be made to fit data from 
SNe la, CMB and BAG [l2l[T3l[Tl[l5l[l6l[I71[T8l[l9]. Further constraints come from 
considering spectral distortions of the CMB [20] and the kinematic Sunyaev-Zeldovich 
effect [21]. 

In the most general scenario, the observer can be located anywhere inside the 
void, in which case observers living off-center will see an anisotropic universe. There 
are several ways to test the level of anisotropy in the data to put constraints on the 
position of an off-center observer. So far, the strongest constraint comes from the 
dipole anisotropy of the CMB. The observed CMB dipole restricts a static observer to 
be located within a couple of percent of the void size [22]. The caveat is - although 
it requires a certain amount of fine-tuning - that a peculiar velocity of the observer 
directed towards the center can cancel the effect of the off-center position so that large 
displacements are allowed. For SNe la, off-center observers will see an anisotropic 
relation between the luminosity distance and the redshift. The SN la data provide 
independent constraints on the observer position that complement those imposed by 
the CMB [23]. Another interesting possibility is to look for cosmic parallax, i.e., a time 
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variation in the angular separation between distant sources induced by the anisotropic 
expansion. Planned space-based astrometric missions (such as Gaia) could measure the 
positions of a huge number of quasars over time with outstanding accuracy and thereby 
establish constraints on the observer position that rival those set by the CMB [23]. 

This paper deals with the contraints on the observer position coming from SNe la. 
The problem has been addressed previously by Alnes & Amarzguioui p3] who found 
that it is possible to obtain a better fit to the data for an off-center observer. The 
constraint on the observer position, however, was found to be fairly weak, partly due to 
the low number of SNe la in the data set employed. Our analysis is similar to theirs, but 
with a few notable differences. In addition to using new and larger SN la data sets as 
well as a different LTB model, we also investigate whether the model can fit the SNe la 
for an off-center observer while simultaneously accommodate the observed CMB dipole. 
We find that it is only possible to obtain a good fit to the data as long as the observer 
is located very close to the void center. 

We would like to emphasise that our analysis is solely aimed at probing the impact 
on the SNe la observations when moving the observer away from the center of the void. 
We employ a simple but plausible LTB model with relatively few parameters to perform 
this test. We do not claim, nor believe, that the LTB model investigated here is able to 
accommodate all the observations coming from other cosmological probes as well. 

This paper is organised as follows. In section Ej we present the basics of the LTB 
model and introduce the specific model employed in this study. In section |3l we introduce 
the SNe la data sets that we use for the analysis. Section H] deals with the case of an 
on-center observer. In section [5l we present the differential equations that govern the 
path and redshift of the photons seen by an off-center observer. The constraint placed 
by SNe la on the observer position is investigated in section |6l In section [7] we look 
at this constraint when the CMB dipole is also taken into consideration. The paper is 
concluded in section [HI 

2. The Lemaitre-Tolman-Bondi model 

A spherically symmetric void can be described mathematically using the Lemaitre- 
Tolman-Bondi metric. 



where the scale function A{r,t) depends on both time and the radial coordinate, and 
k{r) is associated with the spatial curvature. We use primes and dots to denote partial 
derivatives with respect to space and time, respectively. 

In LTB models, the Hubble expansion rate H{r, t) for a matter dominated universe 
can be written as [12] 



A'\r, t) 
1 - k{r) 



dr^ + (r, t) {d9^ + sin^ Odcp^) , 



(1) 




(2) 
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where the density parameters are related by fluif^) + ^k('^) = 1, and the present-day 
values v4o(r) = y4(r, tg) and Ho{r) = H{r,to). Compared to the ordinary Friedmann 
equation of the homogeneous Friedmann- Robert son- Walker universe, all quantities in 
equation ([2]) depend on the radial coordinate. The two arbitrary functions HQ^r) and 
^u{r) define the LTB models and are as boundary conditions independent, allowing for 
inhomogeneities both in the expansion rate and in the matter density. Equation ([2]) can 
be solved analytically using an additional parameter r] [15], 



■^''••''- 211-t.(r)] l'°-''<''>-^-^°''->- 

^°'-"- 2|l-W/- '^""'"''-'''- 

2.1. Gaussian LTB model 

We consider an LTB model in which the matter density parameter Qui^) takes the form 
of a Gaussian density fluctuation, 

nuir) = Qout + (l^in - l^out)e-(^/'^=)' . (5) 

The model has three free parameters, where is the matter density at the center of 
the void, Qout is the asymptotic value of the matter density and is the scale radius 
of the density fluctuation. This density profile is, of course, just a toy model, but it 
provides the basic requirement of a smooth transition from the local to the distant 
matter density, without introducing too many new parameters. 

Our LTB model is more constrained than the general case in the sense that we 
impose that the Big Bang occured simultaneously throughout space by implementing a 
particular choice of Ho{r), 



(6) 



so that the time since the Big Bang is Ibb = f-^^cT^ ^r all observers irrespective of their 
position in space. The model is thus completely specified by only one free function, 
the matter density Quir). Note that the functional form of ^^uir) corresponds to a 
local underdensity also in the physical matter density. The density contrast, i.e., the 
ratio between the local and the asymptotic physical matter density, decreases with 
time, so that the void grows deeper as time progresses. We have introduced the 
pre-factor of 3/2 in equation ([6]) to normalize the age of the universe to that of the 
Einstein-de Sitter universe. The constant Hq sets the scale of the expansion rate and 
determines the age of the universe for the model. Note that Hq does not represent 
the local expansion rate in this model. We marginalize over Hq in the SN la fit, so 
its value is arbitrary in the procedure (as long as rg is in units of Gpc h^^). However, 
for the purpose of establishing an absolute distance scale for the analysis, we choose 
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Figure 1. Distribution of SNe la in the sky in galactic coordinates. Left panel: 
288 SNe la in the SDSS-II data set. SNe la with z < 0.1 are marked with pluses, 
0.1 < z < 0.4 with diamonds and z > 0.4 with squares. Right panel: 397 SNe la in 
the Constitution data set, using the same symbols. 

Hq = 50 km Mpc~iil- This value gives an age of the universe Ibb ~ 13 Gyr and a 
local expansion rate H^lr = 0) ~ 65 km s~^ Mpc"^ for the model. 

Based on predictions from inflation, it is common practice to impose the prior of 
asymptotic flatness, f2out = 1- We will for the most part of the analysis allow f2out to be 
any positive value < 1, i.e., the models need not be asymptotically fiat. In principle we 
also allow for solutions with a local ouerdensity. In this way we can better probe what 
kind of void the SNe la data prefer and also leave open the possibility for curvature in 
the distant universe. We will discuss and compare how imposing asymptotic flatness 
affects the results. 

3. Data sets 

For the analysis, we employ two recent data sets from the literature. The first set is the 
first-year Sloan Digital Sky Survey-II (SDSS-II) Supernova Search data set presented in 
Kessler et al [25]. The compiled data set is based on 103 SNe la at intermediate redshifts 
discovered by the SDSS-II, but also includes previously reported SNe la at low and high 
redshift from other surveys (nearby [26], ESSENCE [27], SNLS [28] and HST [29]), 
bringing the total number of SNe la in the data set to 288. The distance moduli and 
uncertainties were obtained using the MLCS2k2 light-curve fitter. The uncertainties 
include both the observational and the intrinsic magnitude scatter. Figured] (left panel) 
shows the sky distribution of the SNe la in the SDSS-II data set. SNe la with z < 0.1 
are marked with pluses, 0.1 < z < 0.4 with diamonds and z > 0.4 with squares. 
The data set has a comparatively strong emphasis on intermediate redshifts with 132 
SNe la in this range. The SNe la discovered by the SDSS-II cover the redshift interval 



J Since distances depend linearly on Hq, it is easy to convert the quoted values to any other scale 
determined by a different Hq value. 
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z = [0.045, 0.421] and lie on a thin stripe along the equatorial plane. Whereas the 
41 nearby SNe la are fairly evenly scattered across the sky, the 115 high z SNe la 
are confined to small patches. While this sky distribution is not ideal for investigating 
anisotropies in the Hubble diagram, by filling the previously underexplored intermediate 
redshift range, the SDSS-II data set is very interesting for testing LTB models since void 
sizes of the order of 1 Gpc correspond to redshifts in this range. 

The second data set used in our analysis is the Constitution set presented in Hicken 
et al [30]. The data set consists of 397 SNe la and extends the previously available 
Union set [31] (which includes many SNe la from ESSENCE, SNLS and HST) by 
adding 90 SNe la at low reshifts discovered by the CfA3 [32]. The distance moduli 
and uncertainties were obtained using the SALT light-curve fitter. The uncertainties 
include both the observational and the intrinsic magnitude scatter. Figure [T] (right 
panel) shows the sky distribution of the SNe la in the Constitution data set. The set 
contains 141 SNe la at z < 0.1 (pluses), which are evenly distributed across the sky, and 
200 SNe la at 2 > 0.4 (squares), of which the majority are confined to small patches. 
The Constitution set lacks proper coverage at intermediate redshifts. Only 56 SNe la 
occupy the range Q.l < z < 0.4 (diamonds). 

Performing the analysis using these two different data sets is an interesting 
comparison for a couple of reasons. The difference in the redshift distributions affects the 
type of void that the data sets prefer. Meanwhile, the larger number of SNe la and better 
sky distribution of the Constitution set provide stronger constraints on the observer 
position. Another important factor is that the data sets also differ in the light-curve 
fitter used to obtain the distance moduli and uncertainties. Such systematic differences 
can have a large impact on the results of any cosmological test using SNe la [25] . 

4. LTB model for an on-center observer 

In this section we look at the situation when the observer is located at the center of the 
void and investigate the constraints that the data sets infer on the parameters of the 
LTB model. These parameter values are then used as our starting point in subsequent 
sections when the observer is moved off-center. 

4-1. Calculating the luminosity distance for an on-center observer 

An observer located at the center of the void sees a spherically symmetric universe. 
Incoming light travels along radial null geodesies so that the relation between the redshift 
and the coordinates is given by a pair of differential equations [12] , 



dt 

dz 



A'(r,t) 



(7) 



{l + z)A'{r,t) 



dr 

dz 



(8) 



{l + z)A'{r,t) 
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Equations ([7]) and ([H]) determine t{z) and r{z). Combining these with equations (E]) and 
(jl]), we can calculate the angular diameter distance measured by an on-center observer 
as 



The luminosity distance is related to the angular diameter distance according to 



and the distance modulus, which is the difference between the apparent magnitude m 
and the absolute magnitude M, is given by 



4-2. Results 

The best fit LTB model to the SDSS-II data set for an on-center observer has the best 
fit parameters fijn = 0.16, ^out = 0.29 and = 0.7 Gpc (corresponding to z ^ 0.16) 
with = 229.3. In comparison, the best fit fiat ACDM model has = 231.3. For 
the Constitution set the best fit LTB model has = 0.13, Qout = 1 and r^ = 3.5 Gpc 
(corresponding to z ^ 1.02) with = 461.0. Here the fiat ACDM model gives 
= 465.5. Figure [2] shows the results for the two data sets. The top left panel 
shows the functions flui^) and Ho^r) and illustrates well the big difference between the 
best fit models. Whereas the SDSS-II data are best fit with a moderately large void 
with a small change in the matter density, the Constitution set prefers an asymptotically 
fiat model with a very large void. The other panels show the 68.3, 95 and 99% two 
parameter confidence contours (corresponding to Ax^ = [2.30, 5.99, 9.21]) for the model 
parameters, with the SDSS-II contours presented in colour and the Constitution contours 
overplotted with solid lines. There is a clear tension between the data sets, with no 
overlap of the 68.3% contours. This difference in results is partly a consequence of 
the different redshift distributions, but also an effect of the light-curve fitter use 
we impose the prior of asymptotic flatness, flout = 1, the SDSS-II fit becomes notably 
worse, with = 234.6. The local matter density changes to fijn = 0.23 and the scale 
radius is pushed up to a whopping rg = 4.8 Gpc. For the Constitution set this prior 
does not make a difference, since an asymptotically flat model is already preferred by 
the fit. 

The best fits to the data sets are also illustrated in the Hubble diagrams in Figure |3l 
For comparison, we have included the best fit flat ACDM models, and for the SDSS-II 
case we have also plotted the best fit asymptotically flat LTB model. It is quite striking 
how the data differ between the panels. The difference in the best fit parameter values 
for the data sets can be appreciated fairly intuitively from this figure. For the SDSS- 
II data set (left panel), the new SNe la at intermediate redshifts trace out a bump 

§ We note that we see tension also in the best fit fiat ACDM model, where Jim = 0.40 ± 0.04 (68.3% 
confidence limit) for the SDSS-II set and rim = 0.29 ± 0.03 for the Constitution set. 



dA = A[r{z),t{z)] . 



(9) 



C?L = (1 + zfdA , 



(10) 




(11) 
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Figure 2. The best fit on-center LTB models. The upper left panel shows flmir) and 
Ho{r) as a function of the radial coordinate. The other panels show the 68.3, 95 and 
99% confidence contours (for two parameters) for the model parameters f2in, ^^out and 
Ts- The contours are given in colour for the SDSS-II data set and with solid black lines 
for the Constitution set. 

centered around z ~ 0.25, which largely dictates the size of the void. A smaller void 
would place the bump at too low redshift, while a larger void would move the bump 
to higher redshift and make it less pronounced. At higher redshifts, i.e., in the regime 
where flui^) is essentially constant, the curve follows close to a straight line. The low 
value of Qout = 0.29, ensures that it doesn't fall off too fast and misses the high-z 
data bins. The Constitution set (right panel), on the other hand, contains few SNe la 
at intermediate redshifts, but instead features twice as many SNe la as the SDSS-II 
between z = 0.4 and z = 0.9. These give a very broad bump centered around z ~ 0.6 
in the Hubble diagram. The data is best fit with a very large void together with a high 
value for Qont- A smaller void would place the bump at lower redshift, while a decrease 
in the asymptotic matter density would flatten out the curve. 

We also mention that we have investigated the Gaussian LTB model previously for 
the case of an on-center observer and ranked it against the most commonly considered 
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Figure 3. Hubble diagrams for tire best fit on-center models (solid). Distance modulus 
differences in magnitudes are shown with respect to an empty universe. The data have 
been binned for visualization purposes using nAz = 5, where n is the number of SNe la 
in the bin and Az the redshift range. The data points are located at the mean value 
of the redshifts and the weighted mean of the distance moduli. The redshift error bars 
show the standard deviation of the redshifts in the bin, while the distance modulus 
error bars give the error of the weighted mean of the distance moduli within the bin. 
The best fit flat ACDM models (dashed) are included for comparison. For the SDSS-II 
case, we also show the best fit asymptotically flat (Oout — 1) LTB model (dot-dashed). 



dark energy models [33]. We used the SDSS-II data set but also included constraints 
from the cosmic microwave background and the baryon acoustic oscillations in the 
analysis. While the LTB model could provide a fit that was comparable to those of 
the dark energy models, it was not preferred from a model selection point of view since 
extra parameters are required. 



5. Solving the geodesic equation for off-center observers 

We now wish to calculate the luminosity distance to a SN la as seen by an observer 
located off-center in the void. We define the origin of the coordinate system to be at 
the void center and choose the z-axis in the direction of the off-center observer. The 
coordinates of any space-time point in this frame are the time coordinate t, the radial 
coordinate r, the polar angle 6 and the azimuthal angle 0. This means that the spatial 
coordinates of the observer in this frame are r = robs, = and a degenerate 0. Infalling 
photons hit the observer at time to, polar angle 7 and azimuthal angle ^. We need to 
trace these photons back along their trajectories to the source. The redshift of the SN la 
refers to what the off-center observer measures, such that 2;(robs) = 0. 

The system of differential equations that we need to solve in order to obtain the 



coordinates of the SN la as a function of redshift is derived in Appendix A 

dt _ {1 + z) 
dz q ' 



(12) 
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(13) 

(14) 
(15) 

where p = dr/dX, with A being an affine parameter along the geodesic. The constant 
angular momentum J can be calculated as [22] 

J = A(ro,to)sin7 . (17) 

The initial conditions used for the off-center observer are 



dz qA'^ 

dp 1 
dz q 

under the constraint 



A' A' 2 ij2 
-p^ + 



to — tBB 



To = Tohs 



9o = 



cos 7 



(19) 
(20) 

(21) 



A'iro,to) 

The SN la positions in the sky are in the equatorial coordinate system given by the 
right ascension a and declination 6. We need to relate these angles to the coordinate 
system of the void. The polar angle 7 seen by the off-center observer can be obtained 
as 

cos 7 = - sin (3 sin(90° - 6) cos(a + a) + cos f3 cos(90° - 6) , (22) 

where /3 is the inclination angle between the z-axis and the celestial equator's normal 
and a is the angle between the x-axis and the vernal equinox point. The orientation 
angles (3 and a determine the position of the void center and are free parameters in the 
fit. With these definitions, the void center is located at etc = 360° — a and 6c = 90° — (3 
in the equatorial coordinate system. In moving the observer away from the void center, 
we have thus introduced three new free parameters in the fit; the radial displacement 
Tobs and the orientation angles (3 and a. 

Finally, the expression for calculating the angular diameter distance is modified for 
the case of an off-center observer to [23] . 



4 



sin^ e 
sin^ 7 



A'- 



A2(l 



(23) 



where the partial derivatives are obtained numerically in the solution of the geodesic 
equations. 
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Figure 4. Magnitude dipole induced by moving the observer away from the void 
center in the best fit on-center models. The curves show the diff'erencc in magnitude 
for two SNe la with the same redshift but in opposite directions in the sky. Left panel: 
A void with scale radius Vg = 0.7 Gpc {z w 0.18), preferred by the SDSS-II data 
set. Right panel: A void with scale radius Ts = 3.5 Gpc (z w 1-02), preferred by the 
Constitution data set. 



6. Constraining the observer position with SNe la 

Off-center observers will see an anisotropic relation between tlie luminosity distance 
and the redshift for the SNe la. This means that a standard candle with the same 
redshift but in different directions in the sky will have different observed magnitudes. 
The isotropy of the data can be used to establish constraints on the observer position 
inside the void. In this section, we will investigate how far from the center the observer 
can be located. 

6.1. Maximum anisotropy 

To get a sense for how big the effect of being situated off-center has on the SN la 
observations, we can calculate the maximum anisotropy in the form of the magnitude 
dipole. We take two SNe la with the same redshift but in opposite directions in the 
sky, aligned with the off-center observer through the void center. Figure H] shows the 
magnitude dipole, i.e., the difference in magnitude between the two SNe la for three 
different redshifts as a function of the radial displacement of the observer. We have used 
the best fit on-center models for the SDSS-II (left panel), where = 0.7 Gpc {z 0.16), 
and the Constitution data set (right panel), where Vg = 3.5 Gpc {z ^ 1.02). 

The behaviour of the curves is easily understood. For on-center observers, the 
universe is isotropic and the magnitude dipole vanishes. For observers located very 
far from the center, the magnitude dipole becomes less significant. The curves reach 
a maximum at some displacement, depending on the redshift. Figure H] demonstrates 
that SNe la at different redshifts have different constraining power when determining 
the observer position. For the smaller void preferred by the SDSS-II data, the largest 
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Figure 5. The changes in the values relative to the on-center value as a function 
of the observer position. The stars show the angle-optimized values when the void 
model is kept fixed to the best fit on-center parameter values. The diamonds show 
the void model-optimized values when the angles are kept fixed to the angle-optimized 
values. The scale radius of the void is rg = 0.7 Gpc for the SDSS-II set (left panel) 
and = 3.5 Gpc for the Constitution set (right panel). 



anisotropies are obtained for SNe la at 2; = 0.1. The difference in brightness can be 
larger than 0.1 magnitudes if the observer is located a few hundred Mpc from the center. 
SNe la at larger redshifts are less affected by the off-center position. SNe la at low to 
intermediate redshifts will provide the strongest constraints on the observer position in 
this case. In the much larger void obtained for the Constitution set, SNe la at 2; = 0.4 
display the largest anisotropy. An observer located about 1-2 Gpc from the center would 
see a magnitude dipole of around 0.3 magnitudes for these SNe la. For this void, the 
strongest constraints on the observer position come from SNe la at intermediate to high 
redshift. These conclusions also hold for the SDSS-II fit when we impose a prior of 
asymptotic flatness on the model. 

6.2. Results 

We displace the observer from the center and make new fits to the data sets. Figure [5] 
shows the changes in the values compared to the on-center value, as a function 
of the observer position. In the first case (denoted as fixed model) we take the best 
fit on-center model and displace the observer in different directions, i.e., we scan the 
parameters robs, and a for a fixed void model. The stars show the lowest value 
obtained for each value of robs- Note that the scale radius of the void is very different 
for the two data sets, with rg = 0.7 Gpc for the SDSS-II set and rg = 3.5 Gpc for 
the Constitution set. For the SDSS-II (left panel) we find that the fit improves over 
the on-center case out to about robs = 250 Mpc, which corresponds to 36 % of the 
scale radius. The best fit has Ax^ = —0.5 for robs = 100 Mpc and the void center 
at (oc, ^c) = (200°, —50°) in equatorial coordinates or (/c, be) = (308°, 13°) in galactic 
coordinates. For the Constitution set (right panel) we find that the fit improves over 



Supernovae as seen by off-center observers in a local void 



13 



the on-center case out to about robs = 320 Mpc, which corresponds to 9 % of the scale 
radius. The best fit has Ax^ = —1.6 for robs = 160 Mpc and the void center at 
{ac,Sc) = (350°, —50°) in equatorial coordinates or (/c,&c) = (334°, —61°) in galactic 
coordinates. 

As presented in section HI the best fit to the SDSS-II data for an on-center observer 
in an asymptotically flat model is a huge void with r^ = 4.8 Gpc that has Ax^ = 5.3 
compared to the non-flat model. If we displace the observer in such a model, the best 
flt instead has Ax^ = 0.3 for robs = 600 Mpc (13% of the scale radius). A model 
with an off-center observer in an asymptotically flat void can thus provide a flt that is 
comparable to that of an on-center observer in a non-flat void. Note, however, that this 
requires three additional parameters. 

In the second case (denoted as flxed position) we use the orientation angles that 
gave the lowest values in the previous case and instead optimize the void model, i.e., 
we scan the parameters f2in, flout and for a flxed position. Using this approach it is 
possible to improve the flts further. For the SDSS-II (left panel of Figure |5]), the flt is 
improved all the way out to robs ~ 550 Mpc. As the observer is displaced further and 
further from the center, the best flt parameters change to lower values for f2in and fiout, 
and larger values for r^. For the Constitution set (right panel) the flt is improved out 
to robs ~ 400 Mpc. The changes in and are only small and the model remains 
asymptotically flat. 

We would like to point out that we have not performed a full high-resolution 
parameter scan over all six parameters simultaneously, so we expect that it is possible 
to flnd (slighty) lower values than what we have presented here and we refrain from 
giving precise quantitative limits on how far we can be from the center at a certain level 
of confldence. Regardless, our results show that SN la data allow for off-center observers 
and that the flt in fact can be improved for such a scenario. For the Constitution set 
the observer can be displaced ~15 % of the scale radius from the center and still yield 
an acceptable flt to the data. For the SDSS-II set a tolerable flt can be obtained for 
observers displaced all the way out to around the scale radius. However, given that 
three additional parameters for the observer position have been introduced, without 
providing a substantial improvement of the flt, we cannot claim that the off-center 
model is preferred from a model selection point of view. 

7. Constraining the observer position vi^ith SNe la and the CMB dipole 

Being situated away from the center of the void induces anisotropics in the CMB 
temperature [22]. So far we have disregarded this effect in the analysis and focused 
purely on the SN la data. In this section we will continue to investigate the constraints 
on the observer position coming from the SNe la observations while simultaneously 
accommodate the CMB dipole anisotropy measurement. This can be achieved by 
introducing a peculiar velocity of the observer directed to counterbalance the dipole 
induced by the off-center position. Here we will impose that the void model is 
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asymptotically flat, i^out = 1, in order to be consistent with constraints on the spatial 
curvature derived from measurements of the CMB. 

7.1. CMB dipole and peculiar velocity 

The COBE satellite measured the average CMB temperature across the sky to be 
To = 2.725 K and the amplitude of the temperature fluctuations to AT = 3.353 mK |34j . 
The main contribution to the CMB temperature anisotropies comes from the dipole, 
which is in the direction (/, 6) = (264.26°, 48.22°) in galactic coordinates or {a, 6) = 
(168.05°, —7.06°) in equatorial coordinates. In a homogeneous universe, the CMB dipole 
is attributed to a velocity, v^, of the observer relative to the comoving coordinates, 

^ = ^. (24) 
c To ^ ' 

The COBE measurements imply that = 369 km s~^ (corresponding to the net velocity 

of the Sun relative to the CMB). The measured redshifts of the SNe la have been 

corrected for this velocity in the data sets, so that they are given in a frame at rest with 

respect to the CMB. 

In the LTB scenario, the inhomogeneous expansion can be interpreted to cause a 
velocity, fh, of the off-center observer relative to the void center [22], 



rohs ■ (25) 



c c 

In this picture, it is thus the observer's position in an inhomogeneous universe, and 
not primarily a coordinate velocity, that induces the measured CMB dipole. The CMB 
photons arriving from the far side of the void will travel longer through a region with 
larger expansion rate and will thus be redshifted. This means that the CMB temperature 
measured in the direction of the void center will be lower compared to that measured 
in the opposite direction. The void center is then located at (etc, Sc) = (348.05°, 7.06°) 
and the orientation angles are /3 = 82.94° and a = 11.95°. 

It is possible to introduce a peculiar velocity, v^, of the observer, directed towards 
the void center, as a counterbalance in order to allow for larger off-center displacements 
without violating the CMB dipole measurement. The requirement is thus that it is 
the net effect of the off-center position and the peculiar velocity that gives rise to the 
observed dipole, 

Vh + Vp = Vd ■ (26) 

The introduction of a peculiar velocity, directed along the z-axis, will also affect the 
SN la measurements. For a given SN la, the measured redshift z and luminosity distance 
c^L have to be translated into their cosmological counterparts z and c^l according to [35] 

1 + z = (1 -fpCos7)(l + z) , (27) 

dh = di^iz){l + t;pCOS7)(l - t;pCos7)^ . (28) 

We have neglected the peculiar velocities of the SNe la. In the fitting procedure 
we must thus take the measured redshift and obtain the cosmological redshift using 
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Figure 6. The changes in the values for the fit to the SNe la as a function of the 
observer's position. The stars show the values when the static observer is displaced in 
the direction of the CMB dipole in the best fit on-center LTB model. The diamonds 
show the values when the observer also has a peculiar velocity directed to accommodate 
the observed CMB dipole. The arrows indicate the direction of motion, either away 
from the void center or towards it. The vertical dotted line shows the position where 
the peculiar velocity is zero. The scale radius of the void is = 5.0 Gpc for the 
SDSS-II set (left panel) and rg — 3.7 Gpc for the Constitution set (right panel). 



equation fl27|) . calculate the corresponding cosmological luminosity distance, and finally 
translate this value according to equation ( 128|) in order to compare to the measured 
luminosity distance. 

7.2. Results 

We make new model fits for the case of a static on-center observer in an asymptotically 
flat void, using the measured redshifts instead of the redshifts corrected to the CMB 
frame provided in the data set. The best fit to the SDSS-II set has = 231.5 and the 
best fit parameters are = 0.22 and = 5.0 Gpc. For the Constitution set, the best 
fit has = 471.5, with fijn = 0.13 and rg = 3.7 Gpc. 

Using these best fit on-center models. Figure |6] shows how the values change 
compared to the on-center value as the observer is displaced from the center in the 
direction of the CMB dipole. In the first case (denoted as static) the observer has 
no peculiar velocity and the CMB dipole requirement is disregarded. These points 
only serve as a comparison in the plot. In the second case (denoted as moving) the 
observer has a peculiar velocity that perfectly balances the off-center position so that 
the dipole requirement of equation ( l26|) is fulfilled at all values of robs- The arrows 
indicate the direction of the velocity, either away from the void center or towards 
it. The vertical dotted line marks the position where the peculiar velocity is zero. 
Figure |6] demonstrates the power of combining the SNe la data with the CMB dipole 
requirement when constraining the observer position. While the induced CMB dipole 
can always be balanced with an appropriately directed peculiar velocity, this motion 
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will simultaneously affect the quality of the fit to the SNe la. For the SDSS-II set (left 
panel) the peculiar velocity is zero at robs = 28.7 Mpc. Introducing a peculiar velocity 
leads to a deterioration of the fit compared to the static case. The best fit is obtained 
for Tobs = 25 Mpc, which corresponds to 0.5 % of the scale radius. The observer thus has 
a small peculiar velocity directed away from the void center for this fit. While smaller 
displacements can give comparable fits, the value quickly increases as the observer 
is displaced further from the center. For the Constitution set (right panel) the peculiar 
velocity is zero at robs = 23 Mpc. The peculiar velocity leads to an improvement of 
the fit for smaller displacements compared to the static case. The best fit is obtained 
for robs = 6 Mpc, which is only 0.2 % of the scale radius. The observer has a sizeable 
peculiar velocity directed away from the void center for this fit. The value increases 
quickly as the observer is displaced further from the center. 

We conclude that in order to obtain a good fit to the SN la data and simultaneously 
accommodate the CMB dipole requirement, the observer must be located within ~1 % 
of the scale radius. 

8. Conclusions 

We have considered off-center observers in a large, spherically symmetric local void 
described by the LTB metric and investigated the constraints on the observer position 
placed by SNe la. The analysis was performed using two supernova data sets, the 
first-year SDSS-II data set and the Constitution data set, which differ in the redshift 
distribution and the number of SNe la, as well as the light-curve fitter used in order to 
obtain the peak magnitudes. Models with an on-center observer were able to provide 
good fits which had slightly lower values than those obtained for the flat ACDM 
model, but the best fit voids look very different for the two data sets. Whereas the 
SDSS-II data prefer a moderately large void and a low value of the asymptotic matter 
density, the Constitution set is best fit with an asymptotically fiat model with a very 
large void. By displacing the observer from the center we found that the data indeed 
allow for off-center observers. For the SDSS-II set the fit was improved out to 36 % of the 
scale radius and for the Constitution set out to 9 %, the stronger constraint placed by the 
Constitution set as expected from the larger number of SNe la and the more isotropic sky 
distribution. However, for both data sets the improvement of the fit was only marginal. 
Using SN la data alone, we conclude that the observer can be displaced at least 15 % of 
the void scale radius from the center and still give an acceptable fit to the data. These 
conclusions are in good agreement with those reached in Alnes & Amarzguioui [23]. We 
have also performed the analysis using a different parameterisation of the density profile 
[15j with very similar results. 

In the final part of the analysis we also took into consideration that the off-center 
position induces anisotropies in the CMB temperature. While the requirement that 
the induced CMB dipole must be consistent with the measured value can always be 
accommodated by introducing a balancing peculiar velocity of the observer, such a 
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motion simultaneously affects the SN la observations. Using this combination of the 
CMB dipole measurement and the SNe la data, we were able to determine very strict 
constraints on how far from the center the observer can be located. The best fits were 
obtained for an observer located at 0.5 % of the scale radius for the SDSS-II data set 
and 0.2% of the scale radius for the Constitution data set. In order to still get a good 
fit to the SN la data, the observer must be within about 1 % of the scale radius. 

Our more general conclusion is that within the void model, observers have to 
live very close to the center in order for the model to accomodate the data. Besides 
requiring an uncomfortable amount of fine-tuning, such a scenario also constitutes a 
severe challenge to the Copernican principle that we should not occupy a special place 
in the universe. 

In the future, we expect new and better data from many independent cosmological 
probes to put the LTB models to even greater challenges, including new model 
constraints from large SN la data sets with extensive sky coverage also at higher 
redshifts. LTB models may ultimately prove to not be viable cosmological models, but 
they can at least serve as a specific set of toy models to gauge the influence of large-scale 
matter inhomogeneities on the light propagation. Such effects are indubitably present 
as systematic errors in the observations and it is important to be able to quantify their 
significance in order to further advance precision cosmology. 
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Appendix A. The geodesic equations 



The photon paths are governed by the geodesic equation. Alnes & Amarzguioui [22] 
derived the differential equations for the time coordinate t, radial coordinate r and polar 
angle 6 in the LTB space-time. 
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where the paths are parameterized by the affine parameter A and equation flA.3|) has 
been written as conservation of angular momentum J. Note that due to axial symmetry, 
the photon paths are independent of the azimuth angle 0. Furthermore, the equation 
for the redshift z reads as 
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By introducing u = dt/dX and p = dr/dX, we can bring equations (lA.ip and (lA.2p down 
to first order, and the set of differential equations can then be written as 
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Combining Eqs. ( 1A.5I) and (lA.SIl yields 

U = Mo(l + z) , 



(A.5) 

(A.6) 
(A.7) 
(A.8) 

(A.9) 



where we can choose uq = u{X = 0) = —1. 

We want to solve for the coordinates as a function of redshift, and this can be 
achieved since 

d dz d d 
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We can now write down the system of differential equations that we need to solve to 
obtain t{z)^ r^z), 9{z) and p{z), 
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